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Employing thermal dilepton rates based on medium-modified electromagnetic correlation func- 
tions we show that recent dimuon spectra of the NA60 collaboration in central In-In collisions at 
the CERN-SPS can be understood in terms of radiation from a hot and dense hadronic medium. 
Earlier calculated p-meson spectral functions, as following from hadronic many-body theory, provide 
an accurate description of the data up to dimuon invariant masses of about A'f~0.9 GeV, with good 
sensitivity to details of the predicted p-meson line shape. This, in particular, identifies baryon- 
induced effects as the prevalent ones. We show that a reliable description of the p contribution 
opens the possibility to study further medium effects: at higher masses (Af~0.9-1.5 GeV) 4-pion 
type annihilation is required to account for the experimentally observed excess (possibly augmented 
by effects of "chiral mixing"), while predictions for thermal emission from modified u and (p line 
shapes may be tested in the future. 



PACS numbers: 

Introduction. Electromagnetic probes occupy a spe- 
cial role in the study of strongly interacting matter pro- 
duced in energetic collisions of heavy nuclei: once pro- 
duced, photons and dileptons leave the reaction zone es- 
sentially undistorted, carrying direct information from 
the hot and dense medium to the detectors. While the 
ultimate goal in analyzing pertinent experiments is to in- 
fer signatures of QCD phase transitions (chiral symmetry 
restoration and/or deconfinement), the more imminent 
(and relatively easier) objective is to extract medium 
modifications of the electromagnetic (e.m.) spectral func- 
tion P, |3| ■ Recent data from the NA60 collaboration 
on dimuon invariant-mass spectra in In-In collisions at 
the CERN Super-Proton-Synchrotron (SPS) have raised 
the experimental precision to an unprecedented level, 
posing serious challenges to theoretical models. After 
subtraction of late hadronic decay sources ("cocktail"), 
the excess dimuon radiation exhibits features of a (broad- 
ened) peak around the free p-mass (mp~0.77 GeV) with 
a substantial enhancement at both lower and higher in- 
variant masses. Theoretical predictions 0, 0| based on 
an expanding fireball have confirmed that a broadened 
p-meson, as following from hadronic many-body theory, 
is in line with the NA60 dimouns, while a dropping p- 
mass 0, , characterized by a spectrum centered around 
a mass of 0.4-0.5 GeV, is inconsistent with the data^. 
Both scenarios were consistent with earlier dilepton data 
by the CERES collaboration 

The objective of the present letter is twofold: First, 
we improve on earlier (shape-based) comparisons of our 



^ A reinterprctation of the dropping-mass scenario O^within a 
Hidden Local Symmetry approach for vector mesons ^ has not 
been confronted with dilepton data yet. However, any dropping- 
mass scenario will face the challenge of accounting for the large 
portion of strength seen in the NA60 data around the free p-mass. 



predictions to data based on in-medium p spectral func- 
tions. With a slight modification of the expanding fire- 
ball model, we quantitatively reproduce (shape and ab- 
solute magnitude of) the low-mass NA60 data in cen- 
tral In-In collisions. The spectral shape turns out to 
be sensitive to properties of the p in some detail. Sec- 
ond, having determined in this way the contributions 
of the 2-7r piece to the in-medium e.m. correlator, we 
investigate remaining enhancements in the p^ p.~ spec- 
trum. Above M^^~0.9 GeV, using empirical fits to the 
isovector-vector iV) and -axialvector [A) spectra in vac- 
uum, a calculation with the free emission rate underes- 
timates the data. When including medium effects due 
to chiral V-A mixing with an assumed critical tem- 
perature of Tc=175 MeV (consistent with the fireball), 
the description improves. While more precise conclusions 
have to await more accurate data and calculations, our 
estimates indicate that dilepton radiation from heavy- 
ion collisions at the SPS emanates from matter close to 
the expected QCD phase boundary. Finally, we address 
medium effects on the narrow vector mesons ta and </>. 

Thermal Dileptons and Medium Effects. The differen- 
tial rate for thermal lepton-pair production per unit 4- 
volume and 4- momentum can be expressed as Il2l ll^ 



dNu 



a2 L(M2) 



Imn,V^(M,g) /^(qo;r) (1) 



d\.d\ 37r3 M2 
in terms of the retarded e.m. current-current correlator 

n^f^(Q) = i / d'^a:e'«-e(x«) ([J^:,„(x), J,^„(0)]) , (2) 



a final-state lepton phase-space factor, L{M), and the 
thermal Bose function, f^{qo;T). a=e'^ /(Att) =1/137 de- 
notes the fine structure constant and AP^qQ—q^ the 
dilepton invariant mass with energy go a-nd 3-momcntum 
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q. The central objective in the following is the evaluation 
of medium modifications of the c.m. spectral function. 

In the low-mass region (LMR, M<1 GcV), the free 
e.m. correlator is saturated by the light vector mesons 
p, Lu and (f). Our main focus is on the p-mcson (which 
dominates the dilcpton yield in the LMR), but we also 
investigate radiation from (and medium effects on) the uj 
and (j). Although their contribution is suppressed relative 
to the p, it is an inevitable consequence of our assumption 
of the formation of a thermalizcd medium. 

Medium modifications of the p-meson are implemented 
using hadronic many-body theory 0], accounting for in- 
teractions in hot hadronic matter via (i) a dressing of 
its pion cloud with baryon-hole excitations and ther- 
mal Bose factors, (ii) direct resonances on surrounding 
mesons (tt, K, p) and baryons (nucleons, A's, hyperons, 
etc.). The effective interaction vertices (coupling con- 
stants and form factors) have been carefully constrained 
by a combination of hadronic and radiative decay branch- 
ings, by photoabsorption on nucleons and nuclei, and 
by irN—^pN scattering. The resulting p spectral func- 
tions in cold nuclear matter comply with constraints from 
QCD sum rules Q; they have successfully been em- 
ployed to dilepton spectra at full SPS energy 0, and 
to predict an even larger enhancement at lower SPS en- 
ergies 0, 0| . When averaged over a typical space-time 
evolution, the in-medium p width at SPS amounts to 
pmed^35g ]y[gY^ almost 3 times the vacuum value. Close 
to the expected QCD phase the p-resonance has melted, 
rp{Tc)^mp. "Rhosobar" excitations {p^BN~^) lead to 
low-mass strength in the p spectral function that cannot 
be represented in Brcit-Wigncr form. 

Medium effects on the uj and (f> have thus far received 
less attention, especially in the context of ultrarelativis- 
tic heavy- ion collisions (URHICs). For the oj we will 
employ the same approach as for the p The pre- 

dicted average uj width in the hadronic phase of URHICs 
is r^°^~100 MeV 01- For the 0, colhsion rates in 
a meson gas amount to a broadening by ~20 MeV at 
T=150 MeV The dressing of the kaon cloud is pre- 
sumably the main effect for (f) modifications in nuclear 
matter, increasing its width by ~25 MeV at saturation 
density go=0.16 fm~'^ 0. Recent data on (j> absorp- 
tion in nuclear photoproduction suggest even larger val- 
ues |2(1 |. Since a comprehensive treatment of the cj) in 
hot and dense matter is not available at present, we will 
consider the width as a parameter. 

Another important ingredient in our analysis are dilep- 
ton production channels beyond the p, to and (j). For the 
free emission rate these can be inferred from the inverse 
process of e^e~ annihilation, or hadronic r decays as 
measured in decays ji^, enabling a decomposition 
of the (isovector part of) e.m. spectral function into 2- 
and 4-pion pieces. While the former are saturated by 
the p-meson, we fit the latter by an appropriate (onset 




FIG. 1: Free isovector-vector and -axialvector spectral func- 
tions as measured in hadronic r-decays l2lll . compared to fits 
to 3- and 4-7r contributions (the 2-tt piece follows from the p). 

of a) continuum cf. Fig.Q] Rather than evaluating 
medium effects in the intermediate-mass region (IMR; 
1 GeV<A/<3 GeV) in a phcnomcnological approach [23| . 
we here employ model-independent predictions based on 
chiral symmetry. To leading order in temperature one 
finds a mixing of the free V and A correlators [loj . 

UvMq) = (1 - e) Ul^:Aiq) + e n^V (<?) , (3) 

where e=2I(T)/fl (/^=93 MeV: pion decay constant) 
is determined by the pion tadpole diagram with a loop 
integral /(T)=/(d3fc/(2^)3) p{uJk]T)/ujk (wfe: pion en- 
ergy). For £=1/2 (at T^'"'^~160 MeV) V and A corre- 
lators degenerate signaling chiral symmetry restoration. 
Eq. holds in the soft pion limit, i.e., when neglecting 
the pion 4- momentum in Ily^^. A more elaborate treat- 
ment 0, |2^ will broaden and reduce the enhancement 
around M~l.l GeV. To estimate the mixing effect on 
the 4-7r contribution to dilepton emission, we evaluate 
the vector correlator, Eq. with a mixing parame- 
ter £ = /(T, p.n.)/2/(Tc, 0), where critical temperature, 
Tc=175 MeV, and pion chemical potentials (p^>0) are 
in accord with the thermal fireball evolution discussed 
below (TO7r=:139.6 MeV). The 2-pion piece has been re- 
moved as it is included via the p. 

Finally, we account for emission from the QGP, even 
though its contribution at SPS energies is small. We em- 
ploy the hard-thermal-loop improved emission rate for qq 
annihilation |2^, which has the conceptually attractive 
feature that it closely coincides with the rate in hadronic 
matter when both are extrapolated to the expected phase 
transition region. This is suggestive for a kind of quark- 
hadron duality and has the additional benefit that 
the emission from the expanding fireball becomes rather 
insensitive to details of how the phase transition is imple- 
mented (such as critical temperature or "latent heat" ) . 

Thermal Fireball Model. To calculate dilepton spectra 
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we need to specify a space-time evolution of A-A col- 
lisions. Based on evidence from hadronic spectra and 
abundances that the produced medium at SPS ener- 
gies reaches equilibrium |27l | , we adopt a thermal fireball 
model focusing on central In-In with collective expansion 
and hadrochemistry as estimated from data in heavier 
systems. We use a cylindrical volume expansion 



2\2 



VFB{t) = (^0 + v,t) n (r^ + 0.5a^<^) 



(4) 



with transverse acceleration a^=0.08c^/fm, longitudinal 
velocity Vz=c, formation time ro=lfm/c and initial trans- 
verse radius rx=:5.15fm. With hadrochcmical freeze- 
out at (/.t5j'™,T'='^™')==(232,175) MeV and a total fire- 
ball entropy of 5=2630 (using a hadron-rcsonance gas 
equation of state), we have d-/Vch/dy~195~A'part- Isen- 
tropic expansion allows to convert the entropy density, 
s(t)=S/V(t), into temperature and baryon density. The 
evolution starts in the QGP (r(To) = 197 MeV), passes 
through a mixed phase at j"chcm_2T^ g^j-^j terminates at 
thermal frcczeout (rfo~120 MeV). In the hadronic phase, 
meson chemical potentials (/^Tr.if.r;) arc introduced to pre- 
serve the observed hadron ratios Q. When applied to 
dilepton production, the largest uncertainty resides in 
the fireball lifetime (being proportional to the dilepton 
yields), controlled by a± and Tfo, which can be better 
constrained once hadronic data for In-In are available. 
We emphasize, however, that all contributions to the 
dilepton spectrum (QGP, p, uj, (p and 4-7r) are tied to 
the same evolution thus fixing their relative weights. 

Systematic Comparison to NA60 Data. Thermal p^ p~ 
spectra for central In-In are computed by convoluting the 
emission rate, Eq. over the fireball evolution. 
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dNu 
d^a;d^g 



M 

z:^A{M,qt,y), 



(5) 

where A denotes the detector acceptance which has been 
carefully tuned to NA60 simulations The fugacity 

factor, 2;^=e"'^''/^, accounts for chemical off-equilibrium 
in the hadronic phase with n=2,3,4 for the p, uj and 4-7r 
contributions, respectively (for the mixing term in Eq. (O 
we adopt e UY k 0). 

Earlier comparisons [3j of NA60 data to theoretical 
predictions |^ have focused on the contribution from the 
p-meson which dominates in the LMR. While the shape 
of the in-medium p spectral function describes the ex- 
perimental spectra well, the absolute yields were over- 
estimated by ^^20%. This discrepancy can be resolved 
by an increase of the transverse fireball expansion {a±) 
which, on the one hand, reduces the fireball lifetime by 
30% (from 10 to 7 fm/c) and, on the other hand, gener- 
ates harder transverse-momentum (oj ) spectra, which is 
also welcome by preliminary data Consequently, the 
p contribution is reduced, and once QGP emission and 
correlated charm-decays {D,D p^,p~X) are added. 



the spectra in the LMR are very well described, cf. Fig.|21 
Despite the strong p broadening, the p-|-QGP-f charm 
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FIG. 2: NA60 excess dimuons Q in central In-In collisions at 
SPS compared to thermal dimuon radiation using in-medium 
e.m. rates. The individual contributions arise from in-medium 
p-mesons (dash-dotted red line) , 4-7r annihilation with chi- 
ral V-A mixing (dashed blue line), QGP (dotted orange line) 
and correlated open charm (dash-dotted green line); the up- 
per dashed line is the sum of the above, while the solid purple 
line additionally includes in-medium uj and (j) decays 
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sources are insufficient to account for the enhancement 
at M>0.9 GeV. This is not surprising, as 4-7r contribu- 
tions are expected to take over (augmented by a pion 
fugacity factor, e^'^^^'^). Adding the 4-7r piece with chi- 
ral mixing, Eq. nicely accounts for the missing yield 
in the IMR, leading to a satisfactory overall description 
(upper dashed line in Fig. EJ. Going one step further 
still, we argue that in-medium decays of the narrow vec- 
tor mesons, lj and should be included. Whereas their 
decays after freezeout are subtracted as part of the "cock- 
tail" assuming a vacuum line shape, contributions whose 
width goes beyond the experimental mass resolution will 
survive in the excess spectrum. With the predicted in- 
medium Lu spectral function [T^ the agreement between 
theory and data for M=0.7-0.8 GeV seems to improve. 
The (j) contribution is implemented with a two-kaon fu- 
gacity, e^t^Kf^^ and a strangeness suppression factor 7^ 
(7^-0.75 at SPS HI). The </> yield following from the 
fireball model is not incompatible with data, but conclu- 
sions on the spectral shape cannot be drawn at present. 

To better appreciate the relevance of the in-medium 
effects. Fig. 01 shows two scenarios where (i) baryon- 
induced interactions are neglected, and (ii) the vacuum 
e.m. rate is employed. The latter is ruled out; the meson- 
gas scenario, which differs from the full result in Fig. |21 
by factors of 1.5-2 in the LMR, is not favored by ex- 
periment either. Definite conclusions on the in-medium 
p, especially on the baryon-driven enhancement close to 
the dimuon threshold, require a reduction in systematic 
uncertainty, indicated by the open and filled data points. 
Finally, Fig. ^illustrates the signature of chiral V-A mix- 
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M [GeV] 

FIG. 3: NA60 data Q compared to thermal dimuon spectra 
using (i) in-medium p-, to- and (ji-mesons without baryon ef- 
fects (-f QGP+charm-|-in-med.-4-7r; solid purple line), and (ii) 
free p (-f QGP-f charm-f free 4-7r; solid green line). 



4-7r type annihilation) is fixed. Our results confirm the 
prevalent role of a strongly broadened p-meson in the 
LMR, but also suggest substantial medium efi^ects on the 
line shapes of uj and 0. In addition, the IMR might bear 
footprints of chiral Y-A mixing, which would support the 
notion that the matter produced in A-A collisions at the 
SPS is close to chiral restoration. A more quantitative 
connection to chiral order parameters, e.g., by evaluating 
in- medium chiral sum rules [i^l, should be pursued with 
high priority. 
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FIG. 4: NA60 data in the IMR compared to thermal 
dimuon spectra with different implementations of the 4-7r con- 
tribution, using either its vacuum form (lower dashed line) or 
including chiral mixing (upper dashed line), and correspond- 
ing total spectra (lower and upper solid line, respectively). 

ing in the IMR, where the full result of Fig.|5]is compared 
to a calculation with the 4-7r piece in the e.m. rate re- 
placed by its vacuum form (Eq. with e=0). While 
the latter is not incompatible with the data, an identi- 
fication of medium effects requires better precision than 
currently available (both from theory and experiment). 

Conclusion. We have conducted a quantitative in- 
vestigation of the dimuon excess measured by NA60 in 
In(158AGeV)-In. Focusing on central collisions, where 
the notion of thermal radiation is most adequate, we have 
shown that a medium modified e.m. spectral function 
properly accounts for absolute yields and spectral shape 
of the data. While the overall normalization of the spec- 
trum is subject to uncertainties in the underlying fireball 
model (especially its lifetime), the relative strength of the 
different components in the spectrum (p, uj and decays. 
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